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Abstract: In this paper we study finite time future singularities in anisotropic Bianchi type I 
models. It is shown that there exist future singularities similar to Big Rip ones (which appear in 
the framework of phantom Friedmann-Robertson-Walker cosmologies). Specifically, in an ellipsoidal 
anisotropic scenario or in a fully anisotropic scenario, the three directional and average scale factors 
may diverge at a finite future time, together with energy densities and anisotropic pressures. We 
call these singularities “Anisotropic Big Rip Singularities”. We show that there also exist Bianchi 
type I models filled with matter, where one or two directional scale factors may diverge. Another 
type of future anisotropic singularities is shown to be present in vacuum cosmologies, i.e. Kasner 
spacetimes. These singularities are induced by the shear scalar, which also blows up at a finite time. 
We call such a singularity “Vacuum Rip”. In this case one directional scale factor blows up, while 
the other two and average scale factors tend to zero. 

PACS numbers: 98.80.Cq, 04.30.Nk, 98.70.Vc 


I. INTRODUCTION 

The astrophysical observations [l| give evidence that 
our Universe is currently in accelerated expansion. In the 
context of Einstein General Relativity this acceleration is 
driven by an unknown fluid called dark energy 0 , usually 
described by a state parameter w = p/p, with w < —1/3. 
This corresponds to quintessence matter which violates 
the strong energy condition, and the range w < —1 to 
phantom matter, which violates the strong and dominant 
energy conditions. In this latter case we could have the 
scale factor, p and p going to infinity at a finite cosmic 
time in the future. This type of singularity is dubbed Big 
Rip i. This possibility is allowed for isotropic and ho¬ 
mogeneous Friedmann-Robertson-Walker (FRW) models 
by current observational data Q- 

In Big Rip scenarios the curvature invariants R^, 
R^vR^''^ RfivpaR^'^^'^ diverge in the same way as occur in 
the Big-Bang and Big-Crunch singularities H, Q. How¬ 
ever, in the framework of FRW cosmologies there are 
different sorts of finite time future singularities. Accord¬ 
ing to Ref. [ 3 , [ 1 ] the future singularities can be classified 
in the following types: 
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Type / (“Big Rip”) : For t —>• tg, a —>■ oo, p —>■ oo and 
IpI —>■ oo. 

Type II (“Sudden”) : For t ^ tg, a ^ as, p ^ Ps and 
\p\ —>■ oo. 

Type III (“Big Freeze”) : For t ^ tg, a ^ as, p ^ oo 
and IpI —>■ oo. 

Type IV (“Generalized sudden”): For t ^ tg, a ^ Ug, 
p —>• 0, |p| —>■ 0 and higher derivatives of H diverge. 

Type V ( “w-singularities” ): For t — >• G, a — >■ Og, p — >■ 0, 
p —> 0, w ^ oo and higher derivatives of H are regular. 

The quantities tg, ag, ps and pg are constants. 

The tjme II singularity has been studied by several 
authors [9l-[l^ and includes the subcases of the Big Brake 
and Big Boost [l^. In Ref. 0 it was shown that for this 
type of singularity the universe can be extended after the 
singular event. The type III, IV and V have been studied 
in Refs. 0, @,[3] and Q, respectively. 

There are other types of future singularities that can 
appear at a finite time, even when the strong energy con¬ 
dition is satisfied 0, [ij, [ij, 0 . Other interesting types 
of future singularities, but appearing at an infinite time, 
are Little R iv [l3l . Pseudo-Rip [0 and Little Sibling Rip 
of Big Rip [l9l| . It is noteworthy to mention that an at¬ 
tempt to unify future singular behaviors was made in 
Ref. 0, where the authors introduce the Grand Rip 
and Grand Bang/Crunch singularities. 

It is interesting to note that phantom fields are not the 
only way to generate scenarios with Big Rip. Such future 
singularities may be induced, for instance, by fluids with 
an inhomogeneous equation of state [0 or interacting 
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coupled fluids [l^. From the viewpoint of viscous cos¬ 
mological models, in Ref. [2^ it was shown that the bulk 
viscosity induces a Big Rip singularity, and in Ref. 
it was studied a Little Rip as a purely viscous effect. 
Inhomogeneous and spherically symmetric gravitational 
fields, describing evolving wormholes also may exhibit 
a Big Rip singularity during its evolution [2^. Notice 
that in Ref. |ll| an anisotropic and inhomogeneous cos¬ 
mology of Stephani type was found to possess finite-time 
sudden singularities, and in Ref. [l^ specific examples 
of anisotropic sudden singularities in Bianchi type VIIq 
universes were constructed. 

In this paper we extend the study of future singularities 
by considering anisotropic and homogeneous spacetimes 
more general than flat FRW ones. Specifically, we ana¬ 
lyze Bianchi type I cosmologies, allowing us to show that 
the anisotropy of spacetime, by means of the shear scalar, 
may induce future singularities at a finite time, similar 
to Big Rip ones appearing in the framework of phantom 
FRW cosmologies. In order to make analytical progress 
on this topic we shall use some known exact Bianchi type 
I solutions of the Einstein equations, allowing us to han¬ 
dle exact expressions for directional scale factors ai{t), 
shear scalar cr, energy density p and anisotropic pressures 
Pz- 

Our motivation is based on the fact that several stud¬ 
ies on the plausibility of anisotropy in the accelerated 
expanding universe have been performed in the frame¬ 
work of anisotropic dark energy cosmological models. In 
Ref. [ 2 ^ authors found that, in the framework of Bianchi 
I cosmological models, anisotropy is permitted both in 
the geometry of the universe and in the dark energy equa¬ 
tion of state. Additionally, it is worth to mention that an 
anisotropic dark energy model can potentially solve the 
CMB low-quadrupole problem 

The paper is organized as follows. In Sec. II we write 
the Einstein equations for Bianchi type I spacetimes. In 
Sec. Ill we discuss the Kasner vacuum solution and the 
future singularities which may appear during its evolu¬ 
tion. In Sec. IV we find future singularities in anisotropic 
Bianchi type I models filled with a stiff fluid. In Sec. V 
we obtain exact solutions for ellipsoidal universes char¬ 
acterized by a shear scalar proportional to the expansion 
scalar, and filled with isotropic and anisotropic matter 
sources. We show that these spacetimes may exhibit 
anisotropic rip singularities. In Sec. VI we discuss fu¬ 
ture singularities in fully anisotropic Bianchi I cosmolo¬ 
gies filled with an anisotropic, matter source. In Sec. VII 
we discuss our results. 


II. BIANCHI TYPE I SPACETIMES AND 
EINSTEIN FIELD EQUATIONS 

In this paper we consider models belonging to spatially 
homogeneous and anisotropic Bianchi type I spacetimes 
described by the metric 

— al{t)dx^ — a\[t)dy‘^ — a‘^{t)dz^, ( 1 ) 


where ai{t) are the directional scale factors along the 
X, y, z axes, respectively. 

This type of cosmologies is particularly interesting be¬ 
cause it is the simplest generalization of the homogeneous 
and isotropic flat FRW models. 

The Einstein field equations for this metric may be 
written in the following form [ 2 ^: 


9 

(2) 

— 2 H = k{p + p) + 

(3) 

p + SH{p -1- p) = (T • E, 

(4) 

a + 3Hd = E, 

(5) 

where k = SttG, we will consider k = 

1 from here on. 

The average expansion rate i/, the average pressure p, 
the shear vector cr, and the transverse pressure vector E 

are respectively defined as 


// = i (i/i + i/2 + H 3 ) , 

(6) 

p=^{pi +P 2 +P 3 ) , 

(7) 

a,=H,- H, 

(8) 

T,i=pi- p, 

(9) 


where i = 1,2,3. The average expansion rate may be 
written as id = d/d, where the average scale factor d is 
defined by 


a = (010203)^/^, (10) 

and the directional expansion rates are given by 


H, = ( 11 ) 

at 

From Eqs. (|6])-(l9]) we see that the quantities a and S 
satisfy the constraints 

(Ti + a 2 + 0-3 = 0 , ( 12 ) 

E 1 + E 2 + S 3 = 0, (13) 


respectively. 

Additionally we give the definition of some useful 
anisotropic quantities. The shear tensor aab is defined 
by 


^ab — ^a^(c;d)^6 



where 9 = u?^ is the expansion scalar, hab = gab — UaUb 
the projection tensor (for the signature (-I-, —, —, —)), and 
Ua the four-velocity. From this expression we obtain the 
shear scalar given by = (TQh(T“^. 

For the considered Bianchi type I metric © we have 
that the expansion scalar, non-zero shear tensor compo- 
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nents and the shear scalar are given by 

0 = Hi + iJ2 + i73, 

= “3-^2 + g {Hi + H3) , 

O'! = 3 ^^ 2 ), 

^^ = 1 {H^ + Hi + Hi - H 1 H 2 - H 1 H 3 - H 2 H 3 ) , 
respectively. 


where are integration constants for branches + and 
— respectively, and i = 1,2, 3. 

In order to proceed with the analysis, we shall use the 
initial condition Hi{t = 0) = Hq > 0 for the directional 
scale factor oi. This implies that at t = 0 we are imposing 
an expanding scale factor oi. Then, from Eq. (EU we 
obtain that C = and the scale factor along 

x-direction takes the form 


oi = a 


± 

10 



SHot 

I ^<^10 
'^0 



( 22 ) 


III. FINITE-TIME FUTURE ANISOTROPIC 
SINGULARITIES IN VACUUM KASNER 
SPACETIMES 

In this section we study future singularities of 
anisotropic character by considering Bianchi type I space- 
times without matter, i.e. vacuum solutions for the met¬ 
ric © (or Kasner spacetimes). 

By putting p = 0 and Pi = 0 into Eqs. we obtain 

the following four independent differential equations: 



0 (7^ 

377^ = 

(14) 


di + 3Hai = 0. 

(15) 

From Eq. ((ISD 

we obtain 



CiO 

— _3 , 

a3 

(16) 

then cr^ = (cr^Q 

+ O'io + ^ 3 o)/«® S'lid Eq. ([H]) gives 



/ 1 \ 1/3 

a(t) = Ic ±-'/Oaotj . (17) 

Here aio and C are integration constants, and 

CTo = y^2((t^o +crlg + aio(72o), (18) 

where we have used the relation 


0-30 — —CTio — <720. (19) 

Erom Eqs. ®, (HI and (HU) we may write relation 


^ — I 1 ± \ d 

ai \ (To / a’ 


( 20 ) 


which implies that the directional expansion rates Hi are 
proportional to the average expansion rate H. By using 
Eq. (fT71) we obtain for directional scale factors 


^iO 


C ± —\/ 6 ao t 


' i± 


\/6 a 


( 21 ) 


Then, the metric is given by 


ds^ =dH -\l 


mat 


1 ± 




1 -h 


3Hot 


1 ± 


V^cr, 


«( 


1 ± 


V 60 


dx^ - 


! ( 1 I cr2Q \ 

A -0 ) 


dy^- 


1 -f 


3Hot 


[(it 


(q'io+q' 2 o) ^ 




dz^, (23) 


where the constants am have been absorbed by rescaling 
the spatial coordinates. 

Since we are interested in solutions with finite time 
future singularities one should require that 1± < 0. 

This implies that for the positive branch of the considered 
solution the condition 


(Tio < < 0 

V 6 


(24) 


must be fulfilled by coefficients gio and U 2 q, while for the 
negative branch the condition 


cn 

CTlO > —p > 0 


(25) 


must be required. Therefore, in metric (1^ the scale fac¬ 
tor along x-direction exhibits a future singularity at finite 
value of the cosmic time tyr = i ^ > 0. 

Notice that Eq. (EOl) implies that in this scenario we 
have that Hi >0 and H < 0. If we demand that 
the scale factor along p-direction also becomes singular 
at the finite value tyr we must require additionally that 
1 ± < 0, which implies that CT 20 ^ ^ ^ 

the positive branch, and g 20 > > 0 for the negative 

branch. However, it can be shown that simultaneously it 
is not possible to satisfy the conditions gio < ^^id 
tT 2 o < (or gio > and (T 20 > -^). 

Therefore, only the scale factor oi becomes singular at 
finite value of the cosmic time tyy = — (1 i ) > 
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0 , while the other two scale factors 02 and 03 become zero 
at this time (see Fig. [T|). It is interesting to note that 
for all directional scale factors (EU), the corresponding 
expansion rates Hi diverge at tyr- From Eq. o we 
conclude that the same is valid for the average expansion 
rate H. 

It is worth to mention that if we consider C = 0 in 
Eg. dTTl) . then we have C = 0 in Eq. d^Il) and the direc¬ 
tional scale factor ai diverges at tyr = 0 instead of tyy = 
— ^I ± ^. The condition for the occurrence of 

this singularity is the same as before I ± v^ctiq/cto < 0. 
In this sense, the occurrence of the singularity is inde¬ 
pendent of the value of the constant C: if C = 0 we 
must consider t < 0 for the consistency of Eq. m, i-e. 
C ± '/6aQtl2 > 0, if C ^ 0 the consistency of Eg. dTTl) 
C±\^aot/2 > 0 allows us to consider tyy < 0 or tyr > 0 . 

Notice that by defining 


^ 1 I ^ ^_ VGctio _ 

3 y \/ 2 (o-io +cr|o +'^100-20), 

}. (Ij- _ VGa-20 _ 

3 \ \/ 2 ((t^o Tcrlo +'^100-20), 


(26) 


1 


ao ~ ^20 
V^(( 7 l 0 -I- 0-20) 


3 i ^ V^2 ((t^o + o’lo + o’ioo-2o) J ’ 


for the powers of the scale factors in Eq. ((23ll . we obtain 
that 


gi + 92 + 93 = 1, (27) 

9i + 92 + 93 = (28) 

where 91, 92 and 93 are the Kasner parameters. These 
constraints correspond to the conditions for the well 
known vacuum Kasner solution. 

From the Kasner conditions (IT71) and (1281) we note that 
if we arrange the Kasner parameters in increasing order 
9i < 92 < 93, then they change in the ranges [2^ 

- ^ < 91 < 0, (29) 

0 < 92 < I (30) 

I < 93 < 1. (31) 

From these relations we conclude again that if it is 
present a future singularity only one of the scale fac¬ 
tors may blow up, while the other two tend to zero at a 
finite value of the cosmic time. For the particular case 
of an ellipsoidal vacuum cosmology the following param¬ 
eter values must be required: 91 = 92 = 0,93 = I or 
9i = —1/3, 92 = 93 = 2/3. Therefore, only the latter set 
of parameter values allow us to have a future singularity 
for an ellipsoidal vacuum universe (see Fig. [T]). 

In conclusion, due to the anisotropic character of 
Bianchi type I metrics, in the Kasner vacuum solution all 


three scale factors do not exhibit simultaneously a future 
singularity: just one of the scale factors may exhibit such 
a singularity at tyy, while the other two do not. In this 
case the average scale factor (II3 does not exhibit a singu¬ 
lar behavior, becoming zero at tyr = ^. 

We note that this scenario necessarily corresponds to an 
average contracting universe. However, all directional ex¬ 
pansion rates Hi as well as the average expansion rate H 
diverge at and due to the vacuum character of the 
Kasner solutions, the scalar curvature is always zero. 

Because of the absence of matter content, the discussed 
anisotropic future singularities are not similar to any of 
the finite-time singularities listed in the introduction sec¬ 
tion. We shall call such a singularity a Vacuum Rip. 

It should be emphasized that these vacuum rips are not 
produced by fluids violating the dominant energy condi¬ 
tions (DEC) [13, i.e. p > 0 and —p < p < p, as stated for 
FRW cosmologies filled with a phantom fluid. The Kas¬ 
ner vacuum solution satisfies DEC, and by writing the 
Kasner metric in the form where the shear is explicitly 
included, we have shown that future singularities may be 
induced by the anisotropy of the spacetime, by making a 
suitable choice of the model parameters ctio and tT 20 - 


fl; 



FIG. 1 : The figure shows the qualitative behavior of the scale 
factors ai for fully anisotropic (ai 7^ 02 7^ 03) vacuum Bianchi 
type I solutions with crio and <720 satisfying Eqs. (HH) or (I 25 l). 
In the figure are shown two of the scale factors (02 and as) 
which at different rates of contraction tend to zero at the vac¬ 
uum rip time tyr (dotted and dashed lines), while the third 
one, ai, diverges at this time (solid line). In this case the 
future singularity is of anisotropic Cigar Rip type (see TA¬ 
BLE IJ). 


IV. FINITE-TIME FUTURE ANISOTROPIC 
SINGULARITIES WITH A STIFF FLUID 

To elucidate the role of the shear, in the occurrence of 
future singularities, in the presence of matter fields we 
will consider the “toy model” of fully anisotropic Bianchi 
type I spacetimes o, filled with a stiff fluid, for which 
the condition for the powers of scale factors (E3 is still 
valid. This “toy model” is interesting because it allows 
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fl; 



where aj are integration constants. 

By using the initial condition Hi{t = 0) = Hq > 0 for 
the directional scale factor oi we obtain from Eq. (IMl) 

that C = ^ then the scale factor along 

x-direction takes tlie form 


Ui — 


1 + 


3a/12po 


6 crg 


4±- 


a/12po + 6crg ± 6 (Tio 


Hot 


yj 12po+6o-^ 


(35) 


Hence, the resulting metric may be written as 


FIG. 2: The figure shows the qualitative behavior of the scale 
factors Oi for an ellipsoidal vacuum solution given by aio = 
(J 20 = — 1, or equivalently by gi = §2 = 2/3 and ga = —1/3. In 
this case, two of the scale factors at the same contraction rate 
tend to zero at t^r (dashed line), and the other scale factor 
blows up at this time (solid line). Here we have included the 
qualitative behavior of any of the directional expansion rates 
Hi (dotted line), which diverges at t„r- In this case the future 
singularity is of symmetric Cigar Rip type (see TABLE H]). 


ds^ = dt^ 




2 1 

1 ^ -y/lSpo+So-o 

dx' 

7 J 




2 1 4o-20 


j \/l2po+6<Tg 

dy 

1 ) 




2 ^ 4(criQ + g2Q) 

^ l/l 2 P 0 + 6 <^§ ^^2 


(36) 


us to consider finite time future singularities in a cos¬ 
mology, fulfilling energy conditions (p = p, p > 0). The 
other aspect to be considered is that this cosmological 
model allows us to handle exact analytical expressions 
for studying relevant quantities. 

We use the field equations in the form given by 
Eqs. dSl)-®. Since the pressure is isotropic, in this case 
we have that S = 0, and from Eq. o we have for the 
equation of state p = p that the energy density is given 
by p = po/dP- On the other hand, from Eq. ([5]) we ob¬ 
tain the solution (fTHll . and then Eq. ([2]) implies that the 
average scale factor is given by 


a{t) = ± y3po -b ^ 0-2 t j , (32) 

where C is an integration constant and ag is given by 
Eq. (fTSl) . Note that by making po = 0 we obtain the 
average scale factor (II3 discussed in the previous section. 

From Eq. (| 8 ]) we may write for directional expansion 
rates Hi the following equation: 


and the energy density and the pressure are given by 


p = p = 


36iJoVo 


('\/l 2 po + bo-Q ± 


2 ’ 


(37) 


where 


— i j- ^fTio 

3 i/l2po -b 6(7^ ’ 


(38) 


i.e. the power of the scale factor along x-direction. 

In order to induce a future singularity, and considering 
that Hq > 0, we must require that 7 < 0. This implies 
that 


cio < \J 12pg + 6crg < 0, (39) 

for the positive branch, and 

(Tio > i \J 12po + 6crg > 0, (40) 


^ ^ 1 y^l2po -I- 6ao^ 


3 2 C i \/l 2 Po -l- 6 (Tq^ t 

O'iO 


(33) 


C ± y' 3 Po -I-1 fTo^ t 

which implies that the directional scale factors are given 

by 






for the negative branch. 

By taking into account Eq. m we conclude that 


1 


CTio < ^ ( 0-20 - ■\/9o-io + 3po ) < 0, (41) 


for the positive branch, and 


CTio > ^ ( o' 2 o + v 9 0-20 + 8 po ) >0, (42) 


for the negative branch. 
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Notice that by defining the powers of the scale factors 
in the metric (1551) as 


3 ^ 


^ 3 


® = 3T 


21710 

\/l2po + 6crg’ 

2cr20 

\/l2po + 6cr^ ’ 
2 (o'l 0 + O' 2 o) 
\/l2po + 6crg’ 


(43) 


the parameters qi satisfy the condition (EH), indepen¬ 
dently of the values of ctio, <720 and po- Then, the average 
scale factor takes the form 


a = 1 


3a/12po + 6CT| 


y y'Upo + 6crg ± 6crio 
However, now we have that 


Hot 


(44) 


9i + 92 + 93 1 ■ 


(45) 


By putting po = 0 the parameters qi in Eqs. (H51) become 
the Kasner parameters of relations (l26l) . implying that 
the condition (|28l) is fulfilled for vanishing matter. 

Therefore, in the case of Bianchi type I cosmologies 
filled with stiff matter the singularity appears at trs = 
— As in the vacuum case, from Eqs. (03]) we note 
that we can have only one of the directional scale factors 
blowing up together with the energy density and pres¬ 
sure. It becomes clear from expressions ([551)-(0T|) that 
this rip singularity is induced by the anisotropy of the 
spacetime. 

Note that, if we require that two of powers qi are neg¬ 
ative, as it is allowed by the condition (1571) . then po < 0 , 
implying that the energy density becomes negative, vio¬ 
lating the weak energy condition (see Eigs. [5]and|3]). 

The average scale factor (|44p does not exhibit a sin¬ 
gularity at trs, where it vanishes. This anisotropic rip 
singularity appears at trs only for a contracting average 
scale factor. 


0 / 



FIG. 3: The figure shows the qualitative behavior of the scale 
factors at for fully anisotropic Bianchi type I solutions with 
equation of state p = p and satisfying the condition (|27l) with 
qi < 0 and 52 < 0. Here two of the scale factors evolve at 
different rates of expansion and diverge at G (solid and dotted 
lines), while the third scale factor tends to zero at this time 
(dashed line). Note that in this case p < 0 as shown in FIG. 3] 
In this case the future singularity is of pancake rip type (see 
TABLE HI). 



FIG. 4: The figure shows the qualitative behavior of the en¬ 
ergy density (doted line), scalar curvature (solid line) and the 
average scale factor (dashed line) for fully anisotropic Bianchi 
type I solutions with stiff matter. In this case, the energy den¬ 
sity is negative, and together with the Ricci scalar blow up at 
ts, while the average scalar factor tends to zero. 


V. FUTURE ANISOTROPIC RIP 
SINGULARITIES IN ELLIPSOIDAL UNIVERSES 

In this section we consider the evolution of Bianchi 
type I cosmologies with a matter content characterized 
by isotropic and anisotropic pressure. 

Specifically, we consider particular cases of Bianchi 
type I models described by the condition ai(t) = 02 (t)- 
Thus the line element © takes the form 

= dt^ — a\{t){dx^ + dy"^) — a\{t)dz^ , (46) 

which possesses spatial sections with planar symmetry 
and an axis of symmetry directed along the z-axis. The 
functions of the cosmic time ai{t) and 03 (t) are the di¬ 
rectional scale factors along x, y and z directions respec¬ 
tively. The metric (05]) describes a space that has an 


ellipsoidal rate of expansion at any moment of the cos¬ 
mological time, dubbed also Locally Rotationally Sym¬ 
metric Bianchi 1. 

In this case the Einstein field equations ©-(IS]) may be 
written in the form 


p = 

a? , 

9 “T Z , 

(47) 


ai 

Ol 03 


+ 

hi 03 1 03 \ 

(48) 

\ai 

Cl 03 

03 J 

P3 = - 

(2(i + fi), 

(49) 


V 0.1 

aij 


where pi = P 2 and po are the transversal and longitudinal 
pressures respectively. For the metric (1461) the average 
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scale factor is given by a{t) = {af{t) , and the 

average expansion rate takes the form 



a 


1 

3 



(50) 


In order to handle exact solutions to the metric (IMl) . 
we further make the assumption that the scale factors ai 
and 03 are constrained to be given by 


ai{t) = (51) 

where a is a constant. Thus the metric dUl) takes the 
following form: 

= dt^ — + dy'^) — a\{t)dz‘^. (52) 


This metric becomes isotropic for a = 1. It is interesting 
to note that the metric (1521) is characterized by the con¬ 
dition that expansion scalar 6 = {2 + a)H is proportional 
to the shear scalar 

The measure of anisotropy a jO is constant in a number 
of Bianchi-type spacetimes representing perfect fluid cos¬ 
mologies with barotropic equations of state (see [3l| and 
references therein). Given that these models may allow 
nearly isotropic scenarios, they can be used for studying 
the effects of anisotropy in our universe by confronting 
them with observational data. In our case, this condition 
will allow us to work with analytical ellipsoidal cosmolo¬ 
gies exhibiting anisotropic rip singularities. 


A. Anisotropic rip singularities with isotropic 
pressure: p — p 

Let us suppose that Pi = P 2 = Ps = P- Thus from 
Eqs. dm)-® and dm, the relevant metric function 03 (<) 
is given by 


asit) = Cl (1 -I-C21) , (53) 

where ci and C 2 are integration constants. 

We shall rewrite this scale factor by using for the di¬ 
rectional Hubble parameter the condition 


H3{t = 0) = Ho. (54) 

Thus, the scale factor ((551) takes the form 

asit) =ci{l + {2a + l)Hot)T^ , (55) 

and the metric (1551) is given by 

ds^ = dt^ — (1 + (2q; -I- l)ido t) ^“+1 {dx^ + dy'^) 

- {1 +{2a + dz^, (56) 

where the constant Ci has been absorbed by rescaling the 
spatial coordinates. 


In this case the energy density and pressure result to 
be: 


a{a + 2 )Hq 
(1 -f {2a l)ido t) 


(57) 


which means that the isotropic requirement for the pres¬ 
sure implies that the matter filling the universe is char¬ 
acterized by a stiff equation of state. 

From the expression (1551) we see that a future rip 
singularity appears for a < —1/2 when Hq > 0, or 
—1/2 < a < 0 for idg < 0. On the other hand, in order 
to have a positive energy density we must also require 
that a<—2ora>0, which excludes the case Hg < 0. 
For Hq > 0, the scale factor (155)) . energy density and 
pressure blow up at the finite value of the cosmic time 
Us = — (^ 2 a+i)Ho ’ scale factor oi = of becomes 

zero at this time. In this case the average scale factor is 
given by 


d{t) = {l + {2a + l)Hot)^ , (58) 

and does not exhibit a singularity at Gg. 

In conclusion, for the metric (1521) the requirement of 
isotropic pressure implies that the matter content be¬ 
haves as a stiff fluid. The evolution of this cosmology 
exhibits a future singularity for a < —2 {Hq > 0) at 
Us = — (^ 2 a+i)Ho • functions 03 (t), p{t) and 

p{t) blow up at this time, this singularity is similar to 
the FRW Big Rip one but of anisotropic character, since 
at trs the scale factor along x and y directions becomes 
zero, as well as the average scale factor d. From Fq. (1551) 
we note that this scenario corresponds to a contracting 
universe. 

As in Sec. Ill, the anisotropic future singularities are 
not produced by fluids violating the DFC, since in this 
case p = p- This type of singularities is induced by the 
anisotropy of the spacetime, since if shear vanishes, i.e. 
cr = 0 , then the solution becomes the standard isotropic 
FRW cosmology filled with a stiff fluid, which does not 
exhibit any future singularity at a finite value of the cos¬ 
mic time, and only presents the initial singularity or Big 
Bang. 


B. Big Rip singularities with anisotropic pressure 

Let us suppose that the transversal and longitudinal 
pressures are given by 


Pi = uJip, (59) 

P3 = W3P, (60) 

respectively, where wi and luq are constant state parame¬ 
ters. Thus, from Eqs. (SZl), (051) and m we obtain that 

a3(t) = Cl(l -I- C 2 t) + + + l , (0X) 








fl; 



FIG. 5: The figure shows the qualitative behavior of the scale 
factors Oi (dashed and dotted lines) and the average scale 
factor (solid line) for ellipsoidal universes with anisotropic 
pressure. All them blow up at tbr- 



FIG. 6 : The figure shows the qualitative behavior of the en¬ 
ergy density (solid line) and anisotropic pressures pi = p 2 
and P 3 (dotted and dashed lines respectively). In this case, 
all quantities diverge at thr- 


where ci and C 2 are integration constants. By using the 
initial condition ((5^ the scale factor (IFTl) takes the form 

asit) = Cl ^ , (62) 


where 

a -I- 1 

a^LUi -I- -I- 2aa;i -I- a -I- 1 

In this case the metric takes the form: 

ds^ = dt^ — f 1 H-^ tj {dx^ + dy^) — 


(63) 


(64) 


where the constant ci has been absorbed by rescaling the 
spatial coordinates and we have considered i? 3 (t = 0 ) = 
Hq, and the energy density and the pressure pa take the 
form 


P = 


Qf(a -I- 2)Hl 



(65) 


P3 


1 -h 2auJi — a 
1 + a 


P, 


8 

( 66 ) 


respectively. Eq. (| 66 )) implies that the state parameter 
wa is given by 


wa 


1 -I- 2q;wi — a 
1 + a 


(67) 


It becomes clear that, in order to have a positive energy 
density we must require «<— 2 orQ :>0 and for 7 < 0 
the scale factor ( 1621 ) exhibits a future singularity at trs = 
— . At this value of the cosmic time the energy density 

and pressures also blow up. In this case for a < —2 we 
can have one of the scale factors blowing up at trs, while 
for a > 0 all three scale factor may blow up at trs- It is 
interesting to note that the average scale factor, given by 

7 ( 2 ° + 1 ) 

d{t) = fl + ^t^ , (68) 


also may exhibit a singular behavior at trs for 7 < 0 
and a > 0. These inequalities imply, with the help of 
Eq. (1551) . that wi < — ^ 2 a+S ’ hsiice we have that the 
state parameter wi can not be greater than — •\/3/2 for 
any a > 0. For a < — 2 the power of Eq. (l 68 l) is always 
positive and the average scale factor tends to zero at the 
time trs, while the scale factors ai and go to zero and 
oa blows up at this time. This singularity corresponds to 
an axisymmetric pancake rip defined in TABLE U] Note 
that for a > 0 we have H > 0, and for a < —2 we have 
H <0. 

In conclusion, in the case of ellipsoidal universes we 
may have one, or two, or all three scale factors blowing 
up at trs = latter case, it is crucial to 

require a > 0 and a^wi + + 2aoJi -I- a -I- 1 < 0. This 

condition will be realized by requiring 


where 


a±=2- 


a > 1 for wi = — 1 , 

(69) 



r , 

for - 1 < 


a > a- for wi < — 1 . 

(70) 

■(1 “h 2u;i) zb 4:Uj‘f — 3 

(71) 


Therefore, we can have future singularities of cigar and 
pancake rip types (see TABLE |T|. The cigar singulari¬ 
ties may be of anisotropic (ai —>■ 00,02 —>• 0 ,aa —>■ 0 ) 
or symmetric types (oi —>■ 00,02 = oa 0 ), while 
the pancake singularities may be anisotropic and infi¬ 
nite (oi —>■ 00,02 — >■ 00 and oa —>■ 0 ) or of axisymmetric 
(infinite) type (oi = 02 — t 00, oa — >■ 0 ). 

In Figs. [5] and [S] we show the qualitative behaviors of 
the three scale factors, average scale factor, energy den¬ 
sity and pressures for ellipsoidal universes (IMl) . In this 
example the future singularity is of Big Rip type, and the 
universe rips itself apart in all directions at a finite time, 
with diverging energy density and pressures. 


J 
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Initial Singularities 

Anisotropic Rip Singularities 

Type 

Directional scale factors 

Type 

Directional scale factors 

a 

p, \pi\ 

Axisymmetric 

Point-like 

(2i = 0,2 —^ 0, dz —^ 0 





Anisotropic 

Point-like 

Q-i —y 0, (12 —y 0, uz —y 0 





Symmetric 

Barrel 

ui —>■ const., 02 = as 0 





Anisotropic 

Barrel 

ai —>■ const., 02 —>■ 0, 03 —0 





Symmetric 

Cigar 

Oi —>■ 00, 02 = Os ^ 0 

Symmetric 
Cigar Rip 

Ol —)• 00, 02 = 03 ^ 0 

00 

0 or 00 

Anisotropic 

Cigar 

Ol —>■ 00, 02 —>■ 0, 03 ^ 0 

Anisotropic 
Cigar Rip 

Ol —>■ 00, 02 —>■ 0 , 03 ^ 0 

00 

0 or 00 

Axisymmetric 

Pancake 

Ol —>-0, 02 = 03 —>■ const. 

Axisymmetric 
Pancake Rip 

Ol — 0 , 02 = 03 —>• 00 

00 

CXD 

Anisotropic 

Pancake 

Ol —>- 0 , 02 —>■ 00, 03 ^ 00 

Anisotropic 
Pancake Rip 

Ol —>- 0 , 02—^ 00, 03 ^ 00 

cx: 

CX!) 



Axisymmetric 
Big Rip 

Ol —>■ 00, 02 = 03 —>■ 00 

00 

CXD 



Anisotropic 
Big Rip 

Ol —>■ 00, 02 —>■ 00, 03 —>■ 00 

cx: 

CX!) 


TABLE I: In this table all possible initial singularities for Bianchi type I spacetimes are listed (see 0 ). For comparison, we 
also list anisotropic rip type singularities reached at a finite time and described by solutions discussed in this paper. We include 
the directional scale factors ai, shear scalar cr, energy density p and pressures pi. Notice that for vacuum rip singularities only 
symmetric and anisotropic cigar rip types are possible, and only in these cases p and |p| —>■ 0 . 


VI. BIG RIP IN FULLY ANISOTROPIC 
BIANCHI TYPE I COSMOLOGIES 


It is possible to construct a Bianchi type I generaliza¬ 
tion of the ellipsoidal cosmology, exhibiting a future sin¬ 
gularity, with three different scale factors and barotropic 
anisotropic pressures. For instance, let us choose the 
scale factors in the form 

= + ’ (72) 


In this case the average scale factor is given by 


a = 



(Q + g+l)T 

3 


(77) 


We notice that for a = /3 we recover the ellipsoidal cos¬ 
mology of Subsection VB. 

Now we are interested in studying scenarios with 7 < 0 
and a positive energy density. This means that, we must 
require 


a-\- (3 + aj3 > 0. 


(78) 


where si = 07, S2 = Pj and S3 = 7, a and /3 are con¬ 
stants. In this case, from Einstein equations the energy 
density and pressures are given by 


Simultaneously we require that the power of the average 
scale factor in Eq. dZZl) be negative, i.e. 

(a-h/3-b 1)7 < 0. (79) 


P 


{ct cx13)Hq 

(1 + ^)' 


(73) 


From Eqs. (1751) and (1771) we obtain 


a > —1, 


/3> 


a 

1 -I- a 


(80) 


and Pi = LOip (i = 1, 2, 3), where 


Wi 


UJ2 


W3 


1 + /3-7(1 + /3" + /3) 
'){a + l3 + aj3) 

1 -I- a — 7(1 -|- a -I- a^) 
'y{a + f3 + a(3) 
a + P — 7 ( 0 ^ + (3'^ + af3) 
j(a + f3 + a/3) 


(74) 

(75) 

(76) 


It becomes clear that for 7 < 0 at > 0 a 

future anisotropic singularity is present. The character of 
the singularity depends on the values of the constants a 
and 13. From inequations (1801) we obtain two possibilities: 
three divergent directional scale factors for a > 0, /3 > 0 
and two divergent directional scale factors for a > 0, 
/3 < 0 or —1 < a < 0, /3 > 0. 

For positive a and [3 we have Oi —>■ 00, 02 —>■ 00, 
03 —>■ 00, p —>■ 00 and \pi\ —>■ 00. This type of singularity 
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corresponds to an anisotropic Big Rip in TABLE H] For 
q;> 0, /3<0or—l<a<0, /3>0 the singularity 
corresponds to an anisotropic Pancake Rip in TABLE U 
where we also have p —?► oo and \pi \ —>■ oo. 

On the other hand, from Eqs. (Ell)-(ESI) we conclude 
that for positive a and /3 the state parameters wi, uj 2 , W3 
are always negative. In particular, it is possible to have 
a phantom anisotropic cosmology since all < —1, and 
then Pi/p < — 1. Nevertheless, we can not simultaneously 
have three state parameters in the quintessence range 

— 1 < uJi < —1/3, but we can have one phantom state 
parameter and the other two ones in the quintessence 
range. 

In the case of the anisotropic Pancake Rip, the three 
state parameters can be simultaneously positive or in the 
phantom range but all the state parameters can not be 
simultaneously in the quintessence range. In this scenario 
we have two branches: (a) a > 0, /3 < 0 and 7 < 0, 
where we can have one phantom state parameter and 
two state parameters in the quintessence range; and (b) 

— I<a<0, /3>0 and 7 < 0, where we can have only 
one state parameter in the quintesence range and the 
other parameters could be both positives or one positive 
and one in the phantom range. 


VII. FINAL REMARKS 

We have probed that Bianchi type I cosmologies may 
evolve to finite-time singularities, which have anisotropic 
character, since one, or two or even all three scale factors 
blow up at a finite future time. These anisotropic sin¬ 
gularities have certain similarities with that of Big Rip 
type appearing in the framework of homogeneous and 
isotropic phantom FRW cosmologies, in which scale fac¬ 
tor, energy density and pressure become infinite at finite 
future time. 

By considering specific examples illustrating types of 
future anisotropic singularities occurring in the frame¬ 
work of Einstein Bianchi I cosmologies, we study the be¬ 
havior of directional scale factors ai{t), shear scalar a, en¬ 
ergy density p and anisotropic pressures Pi. We show that 
future singularities may be induced by the anisotropy of 
the spacetime. 

In the case of vacuum solutions, i.e. Kasner cosmolo¬ 
gies, only one of the scale factors may exhibit such a 
singular behavior at a finite cosmic time tyr- Tbe other 
two ones tend to zero at this time. In other words, in 
the direction of the scale factor which blows up at tyy 
it follows a super-accelerated expansion until it hits the 
vacuum rip singularity. We call this type of singularities 
vacuum rip. In Figs. [T] and [2] we show the behavior of the 
scale factors for fully anisotropic and ellipsoidal vacuum 
universes. For vacuum rips we have only singularities of 
anisotropic cigar rip type (oi —?► 00,02 — 0,03 —>■ 0) or 
symmetric cigar rip type (oi —>• 00,02 = 03 —t 0), see 
TABLE H 

On the other hand, for non-vacuum Bianchi type I 


spacetimes the scale factors also may evolve to a finite¬ 
time singularity, following a super-accelerated expansion 
until the universe reaches an anisotropic rip singularity at 
which the directional and average Hubble rates, together 
with the shear scalar, energy density and pressure of the 
universe diverge. We have shown that the anisotropy of 
spacetime, by means of the shear scalar, may induce such 
future singularities. 

Notice that this result depends on the initial condition 
we choose for the selected directional scale factor Ui at a 
time Iq (we have chosen i = 1 and to = 0 ). By using the 
condition Hi{to) = Hq the integration constant of the di¬ 
rectional scale factor oi is fixed, therefore expanding and 
contracting scale factors are included in these anisotropic 
scenarios. In this work we chose an increasing directional 
scale factor in order to get an infinite directional scale 
factor at a finite time tg > to- 

To elucidate the role of the shear, in the occurrence 
of future singularities, in the presence of matter fields 
we have considered the “toy model” of fully anisotropic 
Bianchi type I spacetimes ([T]), filled with a stiff fluid, for 
which the condition for the powers of scale factors EH) 
is still valid. This “toy model” is interesting because it 
allows us to consider cosmology, fulfilling energy condi¬ 
tions {p = p, p > 0). The other aspect to be considered 
is that this cosmological model allows us to handle ex¬ 
act analytical expressions for relevant physical quantities. 
The fulfillment of the weak energy condition implies that 
only one of the three scale factors may exhibit such a 
finite-time singularity. Thus, as for the vacuum rip, we 
can have only singularities of anisotropic cigar rip type 
(oi —>■ 00,02 —>■ 0,03 —>■ 0) or symmetric cigar rip type 
(oi —>■ 00, 02 = 03 —>■ 0 ). Accordingly, for p = p > 0 , the 
behavior of scale factors is similar as shown in Figs. [T] 
andm In Figs. [ 3 ] and S] we show the qualitative behavior 
of the solution with two negative powers qi and 52, which 
allows to have negative energy density and isotropic pres¬ 
sure. 

In the case of ellipsoidal universes filled with matter 
characterized by an anisotropic pressure, we may have 
one, or two, or all three scale factors blowing up at tyg = 
— Therefore, we can have future singularities of cigar 
and pancake rip types (see TABLE |T| . The cigar singu¬ 
larities may be of anisotropic (oi —>■ 00,02 —>■ 0,03 —0) 
or symmetric types (oi —>■ 00,02 = 03 —t 0 ), while 
the pancake singularities may be anisotropic and infinite 
(oi 00, 02 00 and as 0) or of axisymmetric (in¬ 
finite) type (oi = 02 —00,03 0 ). In Figs. [S] and [ 3 ] 

we show the qualitative behavior of the three scale fac¬ 
tors, average scale factor, energy density and pressures 
for ellipsoidal universes (l 64 )) . In this example the future 
singularity is of Big Rip type since the universe rips it¬ 
self apart in all directions at a finite time, with diverging 
energy density and pressures. 

It is worth to mention that for a directional scale factor 
Oi(t), exhibiting a singularity at the finite time tg, we 
have that ai(t) —>■ 0 for t —>■ —00, while if a directional 
scale factor vanishes at this finite time tg, then ai(t) —> 00 
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for f — oo. All suitable initial singularities appearing 
for t —?► —oo are included in the TABLE H] From this, it 
becomes clear that the chosen time t = 0 has no physical 
meaning as initial time and only denotes a particular 
moment of the cosmological time during the evolution of 
the considered cosmological models, for which we have 
that —oo <t<ts. In other words, the energy density p, 
anisotropic pressures pi and specific scale factors blow 
up at a time tg from this particular time t = 0. 

We have constructed also a Bianchi type I generaliza¬ 
tion of the ellipsoidal cosmology discussed above, exhibit¬ 
ing a future singularity, with three different scale factors 
and barotropic anisotropic pressures: ai —?► oo, 02 —>■ 00 , 
03 —?► 00 , p —>■ 00 and Pi —00 at tbr- This type of singu¬ 
larity corresponds to an anisotropic Big Rip in TABLE ID 
Note that in this last case the average scale factor and 
the average expansion rate also blow up at tbr- 

Finally, we have shown that it is possible to classify 
finite time future singularities present in Bianchi type 
I models in terms of the evolution of directional scale 
factors. This is done in an analogous way as it is done for 
Bianchi type I models in the case of initial singularities 
[^ . The studied singularities are present when one, two 
or the three directional scale factors blow up at hnite 


time. We note that an anisotropic Big Rip singularity 
(where all the directional scale factors, the average scale 
factor, the average expansion rate, the density and the 
pressures blow up at finite time) is presented in the case 
of ellipsoidal universes with anisotropic pressures and in 
the case of a fully anisotropic universe with barotropic 
anisotropic pressures. Both scenarios correspond to an 
expanding average scale factor, which could represent our 
universe when a small degree of anisotropy is considered. 
We have investigated this last possibility and the results 
will be published elsewhere. 
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